Spin-1 Heisenberg antiferromagnetic chain with exchange and single-ion anisotropics 
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Using density matrix renormalization group calculations, ground state properties of the spin-1 
Heisenberg chain with exchange and single-ion anisotropies in an external field are studied. Our 
findings confirm and refine recent results by Sengupta and Batista, Physical Review Letters 99, 
217205 (2007), on the same model applying Monte Carlo techniques. In particular, we present 
evidence for two types of biconical (or supersolid) and for two types of spin-flop (or superfiuid) 
structures. Basic features of the quantum phase diagram may be interpreted qualitatively in the 
framework of classical spin models. 
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Recently, low-dimensional quantum anisotropic 
Heisenberg antiferromagnets have been shown to exhibit 
the analogue of the supersolid phase^El, usually denoted 
in magnetism as 'intermediate', 'mixed' or biconical 4 
phase, in which both order parameters of the bordering 
antiferromagnetic and spin-flop phases do not vanish. 
Indeed, already some decades ago, in 1956, Matsubara 
and Matsuda^ pointed out the correspondence between 
quantum lattices and anisotropic Heisenberg models, 
when expressing Bose operators by spin operators. Us- 
ing mean-field theory for calculating ground-state and 
thermal properties, supersolid or biconical structures 
have been observed in the uniaxially anisotropic XXZ 
Heisenberg antiferromagnets with additional single site 
terms due to crystal-field anisotropies or with morc- 
than-nearest neighbor interactions^ (note that the 
mean-field approximation of the quantum models cor- 
responds to that of classical models). Such phases may 
give rise to interesting multicritical behavior, especially, 
to tetracritical points^!. 

In the last few years, biconical structures and phases 
in classical XXZ Heisenberg antiferromagnets with and 
without single-ion anisotropies in two as well as three 
dimensions have been studied using ground state consid- 
erations and Monte Carlo technique d 10 * 11 * 12 -! 

Experimental evidence for biconical phases has been 
accumulated over the year d 13 * 14 * 15 * 16 * 17 *. 

The curre nt search for biconical phases in quantum 
magnets^m seems to be partly motivated by the fact 
that they are analogues to the supersolid phased 1 ^* 19 * 20 ! 
Of course, it is also of much interest to study the impact 
of quantum fluctuations on the phases known to occur 
in classical anisotropic Heisenberg antiferromagnets in a 
magnetic field. In the following Note we shall address 
both aspects. 

Specifically, we shall analyze ground state properties, 
T = 0, of the spin-1 XXZ Heisenberg antiferromagnetic 



chain with a single ion anisotropy in a field B. Us- 
ing quantum Monte Carlo simulations, namely stochas- 
tic series expansions, for chains with periodic boundary 
conditions, Sengupta and Batista showed that its quan- 
tum phase diagram at zero temperature displays a field- 
induced supersolid phase™ The model is described by 
the Hamiltonian 

i 

+ D{S?) 2 -BS?) (1) 

where i denotes the lattice sites. For A > 1 and D > 0, 
the exchange and single-ion terms describe competing, 
uniaxial (along the direction of the field B, B > 0, the 
z-direction) and planar anisotropies. Following the pre- 
vious analysis^, we shall deal with the case D = A/2, 
restricting the analysis to the (A, i?/J)-plane. 

To study the model, we here apply density matrix 
renormalization group (DMRG) techniques, yielding ac- 
curate results at zero temperatur e! 21 * 22 ! . In particular, we 
considered chains with open boundary conditions (allow- 
ing the study of fairly long chains) . To monitor finite-size 
effects, the number of sites L ranged from 15 up to 128. 
Usually, a random state was chosen as initial state. To 
get reliable data, typically up to 500 states during 120 
sweeps were kept, with a total truncation error of 10 -7 
and a total energy variance of about 10~ 5 . 

At given chain length L, exchange anisotropy A, and 
field B/J, the total magnetization M = ^ < > 

i 

follows from minimization of the energy E(M, B) — 
Eq(M) — MB, where Eq(M) is the ground state energy 
at B = obtained from the DMRG calculations. Here 
and in the following, brackets, < ... >, denote quantum 
mechanical expectation values. For fixed total magneti- 
zation various physical quantities of interest were deter- 
mined, including the profile of the z-component of the 



2 



magnetization, m, =< Sf >, longitudinal and transverse 
correlators, < S*Sf +r > and, e.g., < SfSf +r >, as well as 
possible order parameters of the various structures and 
phases^. 
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FIG. 2: Magnetization profiles mi for various total magneti- 
zations M in the SS phase in between the IS1 and IS2 phases 
at A= 7 (L = 63). The lower magnetizations belong to odd 
sites i, the upper ones to even sites. The profiles describe 
ground states in the range 10.1 < B/ J < 10.7 



FIG. 1: Ground state phase diagram of Hamiltonian (1) with 
D = A/2, as has been obtained by Sengupta and BatistaP. 
Jj_ in this figure is denoted by J in the text. 

According to the previous analysis^, there are six dis- 
tinct phases in the (A,B/J)- plane at T = 0, as de- 
picted in Fig. 1. The Haldane phase (HP occurs at 
small fields and anisotropics. Of course, the ferromag- 
netic (or normal fluid in quantum lattices) phase (FP) 
occurs for all anisotropics at sufficiently large fields. Fur- 
thermore, there are two solicP or Ising-like phases: The 
usual antiferromagnetic or Neel phase, the TSl'-phase in 
the notation of Sengupta and Batista, and, at rather large 
anisotropies, the 'IS2'-phase with M sa L/2. Actually, 
the IS2-phase corresponds to the (10)-phase in the Ising- 
limit of Hamiltonian (1), the antiferromagnetic Blumc- 
Capel modePS The other two ground state phases are 
the superfluid 6 or spin-liquicP (SL) phase, usually called 
in magnetism the spin-flop phase, and the supersolicP 
(SS) or biconicaP phase. 

Indeed our analysis at selected values of A, 3 < A < 
7, confirms the phase diagram^, see Fig. 1. Moreover, 
we find evidence for two types of biconical as well as 
two types of spin-liquid structures. The evidence stems, 
especially, from the magnetization profiles, m.;. These 
profiles strongly depend on whether the number of sites, 
L, is odd or even, due to the open boundary conditions. 
In addition, finite-size effects may be important, when 
attempting to identify the ground state structures and 
transition points of the infinite chain. 

As illustrated in Figs. 2 and 3, for odd L, the bicon- 
ical structures in between the IS1 and IS2 phases differ 
remarkably from those in between the IS1 and spin-flop 
phases. The first situation is exemplified in Fig. 2 for 
A = 7. We show magnetization profiles m, in the su- 
persolid phase. Close to the center of the chain, i.e. in 



the 'bulk', finite effects have been found to be weak. In- 
creasing the total magnetization M from to about L/2 
one encounters antiferromagnetic and supersolid ground 
states belonging to increasing fields B/J. In the bulk, 
rrii is observed to stay close to one for odd sites i, while 
it changes to almost zero at the even sites on approach 
to the IS2 phase, M i=s L/2. One may describe that be- 
havior in two ways: The spins on the odd sites always 
point in the field, or z, direction, while they turn, start- 
ing from -z direction, in the IS1 phase more and more 
with larger values of M towards the xy-plane on the even 
sites. Alternatively, one may interpret the behavior in the 
framework of the antiferromagnetic Blume-Capel chain. 
At T = 0, there is a direct transition from the antiferro- 
magnetic to the (10)-phase, fixing D and enlargening the 
field. Along that transition line, there is a high degener- 
acy in configurations, where arbitrary fractions of spins 
in the state '— 1' of the antiferromagnetic configuration 
are replaced by spins in the state '0'. These degenerate 
structures seem to give rise to the biconical phase in the 
strongly anisotropic quantum Hciscnbcrg model. 

Attention may be drawn to the weak, but clearly vis- 
ible modulations in the magnetization profiles in these 
supersolid structures, see Fig. 2. 

Let us now consider the supersolid or biconical struc- 
tures in between the IS1 and SL phases. Results for 
local magnetizations, m,, are depicted in Fig. 3 for 
A = 5. Now, on approach to the SL phase, the magneti- 
zation, in the bulk, at the odd sites becomes much lower 
than one, and the difference between the local magneti- 
zations mi in the center of the chain on odd and even 
sites is getting smaller and smaller. Such a behavior 
may be expecte d from classical anisotropic Heisenberg 
antiferromagnets^El] Indeed, for the classical variant 
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FIG. 3: Magnetization profiles in the SS phase in between the 
IS1 and SL phases for various values of the total magnetiza- 
tion M, at A = 5 (L = 63). The lower and upper values be- 
long to odd and even sites, resp. The profiles describe ground 
states in the range 7.0 < B/J < 7.5. 



of Hamiltonian (1), the z-componcnt of the local mag- 
netization has been observed for biconical structures in 
between the antiferromagnetic and spin-flop phases to 
display a behavior similar to that found and seen hertP. 
The spins on both sublattices turn gradually and inter- 
correlated towards their common spin-flop orientation. 

In contrast to the classical variant of the model, the 
transverse components of the spins show no long-range 
(antiferromagnetic) order in the quantum case. Instead, 
the transverse correlations seem to decay algebraically. 

Of course, it would be interesting to clarify whether 
the two types of supersolid or biconical structures are 
separated by a sharp transition and to locate and char- 
actarize that possible transition. Our preliminary results 
suggest that the SL phase seems to disappear at about 
5 < A < 5.5, see also Fig 1. Note that in this region there 
are strong finite-size effects for the magnetization close 
to the center of the chain, and care is needed to discrim- 
inate SS and SL structures. For instance, at A = 5.0, 
we analyzed chains with up to 127 sites, identifying then 
the spin -flop phase. Actually, the possible transition be- 
tween the two distinct supersolid structures may be ar- 
gued to occur in that part of the phase diagram. Detailed 
investigations would require very long chains, and they 
are beyond the scope of this Note. 

In the spin-liquid phase we observe different magne- 
tization pattern for M < L/2 and for M > L/2, at all 
values of A we studied, i.e. independently of the SL phase 
being separated by the IS2 phase or not, see Fig. 1. In 
Figs. 4 and 5, magnetization profiles, mi, illustrate both 
situations. 

As exemplified in Fig. 4, for M < L/2, the local mag- 
netization displays an extended plateau in the center of 
the chain in the SL phase, similar to the behavior in 



FIG. 4: Magnetization profiles in the SL phase with M < L/2 
(L = 63) at A = 3.5. The fields of the corresponding ground 
states are in the range 5.0 < B/J < 5.8. 



the spin-flop phase for classical XXZ Hciscnbcrg anti- 
ferromagnets without and with additional (competing) 
single-ion anisotropy 11 . In Fig. 4 tiny modulations asso- 
ciated with odd and even sites are seen, which, however, 
are even reduced when considering longer chains, possi- 
bly, vanishing for infinite chains. 

Increasing M, M > L/2, but staying in the SL phase, 
the magnetization pattern changes significantly, as de- 
picted in Fig. 5. One first observes a beat-like modula- 
tion about the mean magnetization, as being well known 
from superimposing two sine waves with slightly different 
wavelengths. The wavelength of the envelope of the beat 
is decreasing when enlargening M. Eventually, the mod- 
ulation about the increasing mean value takes on a sim- 
ple (nearly) sinusoidal form, with the wavelength getting 
larger and the amplitude getting smaller when increas- 
ing the total magnetization M. Of course, at M = L, 
the perfect ferromagnetic profile, with Mj = 1 for all 
sites i, is reached. Going from the beat-like to the sinu- 
soidal profiles, one increases systematically the average 
distance between successive extrema, reflecting, presum- 
ably, a quasi-continuous increase of the winding number 
of the modulation. Thence, it seems tempting to sug- 
gest the SL phase at M > L/2 to be of incommensu- 
rate type, while it seems to be of commensurate type 
for M < L/2. A similar distinction for the spin-flop 
phase may have been observed before in different parts 
of the phase diagram of Hamiltonian (ljpSI It is worth 
mentioning that the amplitude of the modulations shrink 
somewhat for longer chains, and it would be interesting 
to analyze whether such, presumably, Fricdel like oscil- 
lations are still present in the infinite chain. 

As in the case of the biconical phase, the transverse 
components of the magnetization show no long-r ange or- 
der in the SL phase of the quantum mode£UM23 in con- 
trast to the situation in the classical variant of Hamilto- 
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Again, a more detailed analysis, considering carefully 
finite-size effects, is desirable, but well beyond the scope 
of our present investigation. 

To summarize our findings, we have studied the S-l an- 
tifcrromagnetic XXZ Hcisenberg antiferromagnetic chain 
with a competing, planar single-ion anisotropy. Using 
DMRG calculations, we confirmed and refined the ground 
state phase diagram obtained recently by Sengupta and 
Batista applying Monte Carlo techniques. In particular, 
we presented evidence for two distinct types of bicon- 
ical or supersolid structures and for two distinct types 
of spin-flop or superfluid structures, in the finite chains 
we studied. We compared the findings with results on 
related classical magnets. 



FIG. 5: Magnetization profiles ra» in the SL phase with M > 
L/2 (L = 63) at A = 3.5. The fields of the corresponding 
ground states are in the range 6.8 < B/J < 10.7. 

nian (1). 

Spin-flop phases with beat-like or sinusoidal modula- 
tions in m, do not occur in the classical variant of the 
Hamiltonian (1). 
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